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It has recently shown that in supersymmetric left-right models with automatic R-parity con- 
servation, the theory below the Wr scale is given by MSSM with massive neutrinos and a pair of 
doubly charged superfields with masses in the 100 GeV range ( with or without an extra pair of 
heavy Higgs doublets (M> 10 TeV) depending on the model). In this paper we study the unification 
prospects for such theories and their phenomenological implications for collider experiments. We 
study two versions of the theory, one with supersymmetry breaking transmitted via the gauge and 
another where the same occurs via gravitational forces. We point out that looking at multi r final 
' states can considerably constrain the parameter space of the model. 
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Supersymmetric left-right models (SUSYLR) where the SU(2)r gauge symmetry is broken by triplet Higgs fields 
with B — L = 2 have many attractive features: 1) they imply automatic conservation of baryon and lepton number [Q, 
a property which makes the standard model so attractive, but is not shared by the minimal supersymmetric standard 
. model (MSSM); 2) they provide a natural solution to the strong and weak CP problems of the MSSM [^j; 3) They 
O^l ' yield a natural embedding of the see-saw mechanism for small neutrino masses || where the right handed triplet 
field (A c ) that breaks the SU{2)r symmetry also gives heavy mass to the righthanded Majorana neutrino needed for 
implementing the see-saw mechanism. 

In order to cancel anomalies as well as to maintain supersymmetry below the SU(2)r scale (vr), one needs a pair 
of fields A c © A c with B — L = =f2 respectively. An important distinguishing characteristic of these Higgs multiplets 
\ is that they contain doubly charged Higgs bosons and Higgsinos in them that remain as physical fields subsequent to 
symmetry breakdown. It has recently been shown that the vacuum of the theory may or may not conserve R-parity 
O i' [4|. If however we require the ground state to conserve R-parity , one must include higher dimensional operators 
or additional Higgs fields which break parity 0,0. In this case, ]t]|| the model in its simple versions, always 
' predicts that some of the doubly charged fields mentioned above are massless in the absence of the higher dimensional 
operators (HDO). This is independent of whether the hidden sector supersymmetry breaking scale is above or below 
the Wr scale. In the presence of HDO's, they acquire masses of order ~ v R /Mpi. Since the measurement of the 
Z-width at LEP and SLC implies that such particles must have a mass of at least 45 GeV, this puts a lower limit on 
the Wr scale of about 10 10 GeV or so. This result is interesting since Wr masses in this range also lead to neutrino 
masses expected on the basis of current solar and atmospheric neutrino experiments. If we take the Wr mass to be 
close to this lower limit (say of order 10 10 — 10 11 GeV), it implies the masses of the doubly charged particles in the 
100 GeV range. There are two bosonic and two fermionic particles of this type. The rest of the particle spectrum 
below the Wr scale can be same as that of the minimal supersymmetric standard model (MSSM) with a massive 
neutrino or it can have an extra pair of Higgs doublets in the 10 TeV range depending on the structure of the model. 
It is the goal of this paper to explore the constraints on the parameters of the model and suggest tests in the e + e~ 
and pp collider. 

Our main results are that for a large range of parameters in this model the stau (the superpartner of r) is the 
lightest of the sleptons due to the renormalization group running arising from the A ++ t~t~ coupling which is not 
very constrained from phenomenology. (Note that the A ++ couplings to other leptons are severely bounded by the 
recent PSI results on muonium- antimuonium oscillation S). As a result, tau lepton final states, generated from the 
production of the A ±=t and its fermionic part, provide a crucial signature of this class of models. For instance we find 
that detection of final states of type t~t~t + t + 77 plus missing energy or t~t~t + t + with or without missing energy 
in both pp and e + e~ collsion will provide test of these models. Thus non-observation of such signals will significantly 
reduce the domain of allowed parameters of this model. We point out the difference between the multi r signals of 
this model and the same type of signals that appear in the other models e. g. conventional GMSB type. We will 
show that the allowed region of parameter space is larger when we choose GMSB type of theories. We also discuss 
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the unification prospects of this model in different SUSY breaking scenarios. 

This paper is organized as follows: in section 2, we review the arguments leading to the existence of light doubly 
charged fields despite a high Wr scale in the context of a simple model and discuss the low energy interactions of 
these fields; in section 3 we discuss gauge unification in these models; in section 4, We discuss its parameter space 
and experimental signals in gravity mediated SUSY breaking scenario; in section 5, we discuss the same for the gauge 
mediated SUSY breaking scenarios and in section 6 we present our conclusion. 

II. OVERVIEW OF THE MODEL 

In this section, we present a brief review of the arguments leading to the existence of the light doubly charged Higgs 
fields in the SUSYLR model. In order to give our arguments, we start by giving the basic features of the model, which 
is based on the gauge group SU(2)l x SU(2)r x U(1)b-l x SU(3) c . In Table I, we give the particle content of the 
model. We will suppress the SU(3) C indices in what follows. 



Fields 


SU(2) L x S\J(2) R x U(1)b_ l 
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(1,1,0) 
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Table l:Field content of the SUSY LR model; we assume that S 

is odd under parity; U and V denote the SU(2)l,r transformations respectively. 

The superpotential for this theory is given by (we have suppressed the generation index): 

W = h«g T T 2 $ l T 2 Q C + h ; (4) i T T 2 $ 4 T 2 i C 

+ i(iL T T 2 AL + i c L cT T 2 A c L c ) 

+ M A [Tr(AA) + Tr(A c A c )] + AS(AA - A c lY) + fi s S 2 

+ fi ij Tr(T 2 <f>jT 2 $ j ) + W NR (1) 

where Wnr denotes non-renormalizable terms arising from higher scale physics such as grand unified theories or 
Planck scale effects. 

W NR = ^[Tr(A c A c )] 2 /2 + BTr(A c A c )7>(A c A~ c )/2 (2) 

where A and B are of order l/M P i anck . 

We will work in the vacuum which conserves R-parity. The Higgs vevs then have the following pattern: 

<*>=(» °) :<at> =(oS) o 

Similar pattern for < A c > is assumed. 

Using Eq. 1, one can give a group theoretical argument for the existence of light doubly charged particles in the 
supersymmetric limit as follows. Let us first ignore the higher dimensional terms A and B as well as the leptonic 
couplings /. It is then clear that the superpotential has a complexified U(3) symmetry (i.e. a t/(3) symmetry 
whose parameters are taken to be complex) that operates on the A c and A c fields. This is due to the holomorphy 
of the superpotential. After one component of each of the above fields acquires vev the resulting symmetry is the 
complexified U(2). This leaves 10 massless fields. Once we bring in the D-terms and switch on the gauge fields, six 
of these fields become massive as a consequence of the Higgs mechanism of supersymmetric theories. That leaves 
four massless fields in the absence of higher dimensional terms. These are the two complex doubly charged fields. 
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Of the two non-renormalizable terms A and B, only the A-term has the complexified U(3) symmetry. Hence the 
supersymmetric contribution to the doubly charged particles will come only from the B-term. It is then clear that 
the masses of the doubly charged fields are of order v%/Mpi. Requiring that these masses satisfy the Z- width bound 
then implies that vr > 10 10 — 10 11 GeV or so. In this paper we will assume that vr is at the lower bound value so 
that the doubly charged fields are accessibly to the existing collider experiments. Note incidentally that although the 
leptonic couplings do not respect the above mentioned symmetry, they are unimportant in determining the vacuum 
structure as long as R-parity is conserved and hence they do not effect the doubly charged field masses. 

Let us now give an explicit calculation of the masses of the doubly charged fields in the supersymmetric limit using 
the superpotential in Eq.(l). Let us write down the F-terms for the S, A and A c terms: 

F s = 2fi s S + A(AA - A c A~c) 
F A = (AS + M A ) A 
F A a = (-AS + M A )~E? (4) 

If the effective supersymmetry breaking scale is below the Wr scale, then these F terms must vanish. It is then 
clear that if we choose the A c and A c vev's (denoted by Vr and vr) to be nonvanishing (and they are equal in the 
supersymmetric limit), then we must have < S >— Ma/A. This implies that the A and A vev's vanish and the 
masses of these fields are of order 2M A . Thus the left triplet fields decouple from the low energy spectrum. It is then 
easy to see from the superpotential (in the absence of the A and B terms) that all the particles in the superfields A c 
and A c are massless in the limit of exact supersymmetry. One linear combination of the neutral fields and another 
ofsingly charged fields disappear due to the Higgs mechanism. The remaining singly charged and neutral Higgs fields 
pick up mass of order of vr and disappear from the low energy spectrum. 

The story of the doubly charged fields is however very different in this theory as has been shown in Ref. ||. Once 
supersymmetry breaking is turned on but the higher dimensional terms A and B are excluded from the analysis, these 
fields acquire negative mass-squares signalling the breakdown of electric charge. This problem is cured as soon as the 
A and B terms are included. The doubly charged fields (all of them) then acquire masses of order v^/Mpi and the 
vacuum becomes charge conserving. 

Let us now discuss the Higgs doublet spectrum of the model at low energies. At the Wr scale, one generally takes 
two bi-doublet fields <fi's to make the model realistic. In order to get the MSSM at low energies, one must decouple one 
pair of H u and from the low energy spectrum. This has been called doublet-doublet splitting problem in literature. 
In the model without HDO contributions, it is clear from the superpotential in Eq. (1) that doublet Higgsino matrix 
is symmetric: 

M H = ( ^ ^ ) (5) 

If we now do fine tuning to get one pair of H u c i at low energies, the H u d appear as identical combinations of the 
doublets in </>,;'s. As a result, at the MSSM level, we have proportionality of the M u and Md leading to vanishing 
CKM angles. This result holds even if we increase the number of bi-doublets arbitrari ly a nd uses only the fact that 
bilinear mass matrix ynj is symmetric. This in fact raises the interesting possibility |1C[ | that all mixing angles in 



the quark and lepton sector may arise purely out of radiative corrections involving the soft breaking terms 1 1 1 . An 
advantage of this version of the model is that there are no new flavor changing effects other than those from the usual 
supersymmetric sources |T^| . 

On the other hand, one may choose the Hij parameters of this model to be of the order of electroweak scale so that 
the low energy model is not exactly the MSSM but rather the two Higgs pair extension of MSSM. The phenomenology 
of these models are very similar to the previous case except that there are new contributions to the flavor changing 
neutral current effects in this model similar to those in the nonsupersymmetric left-right models [fl3| which puts a 
lower limit on the masses of the second pair of Higgs doublets to be in the 5-10 TeV range. As a result they will 
essentially decouple from the low energy spectrum. 

Our results are independent of which of the above choices for the Higgs sector is made, except that unification 
discussion applies only to the second version. 



III. GAUGE UNIFICATION 



The presence of the doubly charged fields at low energies distinguishes the gauge coupling evolution in this model 
from the MSSM and one might expect that one will lose the unification property. It however turns out that the gauge 
couplings do unify in this model, albeit at a lower scale as we see below, for the case which has two pairs of Higgs 
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doublets at the weak scale. The couplings evolve according to their respective beta functions. As just mentioned, 
below the vr scale, we assume four Higgs doublet fields instead of the usual two of MSSM and we assume that one 
of the doublet pairs has mass of 10 TeV. They lead to a trivial modification of the beta function below vr. The beta 
functions above the vr scale are given below for one loop : 
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and in the following equation for the twodoop: 
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where i = U (1)b-l, SU(2)l 1 SU(2)r, SU(3)c respectively in the matrices, Np is the number of fermion generations. 
Np = 3 always and is the number of bidoublets which we take to be 2. We also take tia (A + A) and tia<= ( A c + A c ) 
to be 1. Since the A ±=t leaks down to the weak scale, we need to include its contribution to the running of the gauge 
couplings in between the weak and the intermediate scale. Since this field A has only hypercharge quantum number 
under the SM representation, the hypercharge gauge coupling RGE gets an extra term of 24/5 in one loop and in the 
two loop, the hypercharge squared elements gets and additional factor 72 x 16/25. 

We see from figure la. that the gauge couplings unify at a scale <~ 10 11 GeV and the intermediate scale is 
~ 10 8 GeV. We take a c = 0.118, a = 1/128.7 and sin 2 9 W = 0.2321 at the weak scale. The low unification scale 
implies that the proton decay constraint would rule out groups like SO(10). However the final unifying group can be 
SU(3) x SU(3) x SU(3) or any group that conserves baryon number. 

In the supergravity motivated models, the gauge unification is necessary in order to have gaugino mass unification. 
The masses of all the sparticles can then be determined in terms of the parameters e.g. universal scalar mass mo at 
the unification, universal gaugino mass mi/ 2 and the trilinicar coefficient in the potential A's. We will assume the 
A's to be 0. We choose to assume the universality of too at the unification scale rather than at the Planck scale for 
two reasons: first is that we do not know the theory above the unification scale and the nature of the evolution of 
the parameters is obviously dependent on those details. The second reason is that the experimental signature we are 
interested in involves only the tau lepton and its partner and not the other super particles of the theory. So even if 
we assumed universality of scalar masses at the Planck scale which would necessarily imply some splitting between 
generations due to running between the Planck scale to the unification scale, our final conclusions will be uneffected 
by this. 

If we use GMSB (gauge mediated supersymmetry breaking models) models where SUSY breaking is communicated 
to the observable sector by gauge mediation, the soft susy breaking scalar and the gaugino masses are generated 
at a scale ~ 10 5 GeV by gauge interactions. In this type of models gauge unification is not necessary in order to 
have unified gaugino mass, since they are generated in one loop. In our calculation for the GMSB cases we do not 
consider any gauge unification. However in figure lb we show one example where we have the gauge unification with 
a messenger sector composed of: one copy of Dr + Dr, Ur + Ur and Ll + L^. The couplings unify at 10 10 5 GeV, 
the intermediate scale is at <~ 10 5 2 GeV and the ctR/otB-L is 1/2 at the intermediate scale. 
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FIG. 1. Gauge coupling unification (two loop) in the Supergravity case is shown 
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FIG. 2. Gauge coupling unification (two loop) in the GMSB case is shown 

IV. GRAVITY MEDIATED SCENARIO 

In the standard gravity mediated scenarios, one starts with a universal mass-square for all scalar components of 
the chiral superfields at the Planck scale and they are then extrapolated to the weak scale to determine their masses 
Q. For those fields with large Yukawa couplings such as the Higgs doublet H u , top squark, the weak scale value 
is significantly lower than the planck scale one. (In fact, for the H u turning negative gives rise to the celebrated 
phenomenon of radiative electroweak symmetry breaking). For other squarks, the gaugino mass contribution has the 
effect of increasing their value over the Planck scale value. On the other hand for the sleptons, the change between 
the weak scale and Planck scale value is not very significant since they neither have large Yukawa couplings nor do 
they have strong interactions. 

Turning now to the SUSYLR model, as mentioned before, we will assume universality at the unification and in 
contrast with the MSSM, the effective theory below the right handed scale, contains the following new coupling in 
the superpotential: 

W extra = f i A c±± kk (8) 

which will have a major impact on our spectrum at the weak scale. 

In writing the above coupling, we have used the fact that experimental limits on lepton flavor changing processes 
such as fi — > 3e and r — > 3e imply that with i ^ j are very small compared to the diagonal couplings fa and have 
therefore taken the liberty to simply drop the off-diagonal couplings. In what follows, we will denote A c = A. This 
interaction gives rise to the process /i + e~ — * /i~e + with a strength G M _^ ~ Jyf^r- Recent PSI experiment j| has 

yielded an upper limit on G M _ A j < 3 x G F x fCT 3 . For M A = 100 GeV, this implies that /i/ 2 < 1.2 x 10~ 3 . Thus we 
expect each of the couplings to be less than . 1 barring pathological situations where only one of the couplings bears 
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the brunt of the constraint. On the other hand there is no such constraint on from experiments. So we will choose 
it to be of order 0.5. 

The first implication of the relatively large f% is on the weak scale value for the mass of the A-boson. As shown 
in ||, at the vr scale, both the A and its fermionic partner A (A-ino) have nearly the same mass. The bosonic 
component however runs faster than the fermionic one. As a result, at low energies we can expect that Ma < M A . 
This has important implications for phenomenology 

A second consequence of the possible large fz is that the t c mass is drawn down to smaller values at lower scales. 
In order to calculate the mass spectrum we use the fermion masses, mixing angles, /'s and the magnitude of the gauge 
couplings as inputs at the weak scale. We use m t = 175 GeV and mi, = 4.5GeV. We then use the RGE's shown in 
the appendix and in the reference to determine the Yukawa and the gauge couplings at the GUT scale. Then we 
use the universal boundary conditions and use the RGE's for the soft SUSY breaking masses and the couplings to 
determine the masses at the weak scale. We determine the parameter /i from the condition that electroweak symmetry 
is broken radiatively. We choose the sign of /i to be negative, since the other choice will give large b sj rate. The 
lighter stau mass will be much smaller than the other sleptons (even when tan/3 is small) due to the presence of the 
new coupling f%. In figure 3 we show the mass contours of the lighter stau (dotted line), the A ±± (solid line) and the 
lighter selectron(dot-dashed line) in the too and tan/3 plane for m 1 / 2 — 180GeU. We choose A ±=t mass ~ 100 GeV at 
the Weak scale (90 GeV at the intermediate scale). As mo increases the A ±=t mass decreases due to the subtractive 
effect in the RGE originating from the larger soft SUSY breaking scalar mass. We can see from the figure that the 
bound on A ±=t mass rules out the upper range of mo. The lower bound on lighter stau (fi) mass on the other hand 
rules out lower range of too . 

The latest bound on lighter stau mass is about 57 GeV jl6) . In this scenario A mass is too high to be pair produced 
at LEPII. We also observe (due to the smallness of /i,2 compared to /a) that the lighter selectron mass is much higher 
than the fj mass even when the tan/3 is small for the same Too and mi/ 2 values. This will differentiate between the 
final states of chargino pair production in MSSM and in this model. In fig 4 we show the same mass contours for 
A ±± ~ 80 GeV at the weak scale (70 GeV at the intermediate scale). Lower A ±=t mass indicates lesser effect on the 
fx mass from the new interactions and consequently more parameter space for lower Too, however low A mass rule 
out more parameter space from the higher too range. In this scenario the As can be pair produced at LEPII. In Fig 5 
we exhibit a scenario where A mass is ~ 120 GeV at the weak scale (110 GeV at the intermediate scale). The lightest 
neutralino (xi)in all these scenarios are ~ 49 — 57 GeV and the lightest chargino mass is around 80-90 GeV. If we vary 
the 7774/2, e.g TO1/2 = 150GeV, the x? becomes 33-43 GeV and the lightest chargino becomes 57-67 GeV. the scalar 
masses also get reduced by 10 -20 GeV. In Fig 6 we show a scenario where the third generation coupling is smaller 
/3 ~ 0.25. As one can expect, the effect of the new couplings in the stau mass is reduced. For low tan/3, as in the 
conventional SUGRA mode, the t\ and the lighter selectron mass (en) become very close. The A boson mass lower 
than 88 GeV appears for to )270 GeV, where as in Fig. 3 the same mass contour appears at Too) 140 GeV (In both 
the cases the m^/i at the GUT scale and M^ at the intermediate scale are same. From the analysis of the parameter 
space we surmise that most of the allowed parameter space for large fa can be searched in the present colliders. 




100 150 200 250 300 



m„(GeV) 

FIG. 3. The mass contours for M^±± — 90GeV at the intermediate scale and fz — 0.5 are shown. The fi mass contours 
(dashed line) from left to right depict 45, 60, 80 and 100 GeV masses, the A ±± mass contours (solid lines) from left to right 
depict 88,70 and 50 GeV masses, the en (in and fin masses are same) mass contours (dot dashed line) from left to right depict 
100 and 140 GeV masses 
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FIG. 4. The mass contours for M^±± — 70GeV at the intermediate scale and fz = 0.5 are shown. The fi mass contours 
(dashed line) from left to right depict 45, 60, 80 and 100 GeV masses, the A mass contours (solid lines) from left to right 
depict 80, 65 and 50 GeV masses, the in mass contour (dot dashed line) depicts 100 GeV mass 
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FIG. 5. The mass contours for M^±± = HOGeV at the intermediate scale and fs = 0.5 are shown. The fi mass contours 
(dashed line) from left to right depict 45, 60, 80 and 100 GeV masses, the A ±± mass contours (solid lines) from left to right 
depict 108, 95 and 70 GeV masses, the Br mass contours (dot dashed line) from left to right depict 100 and 140 GeV masses 
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FIG. 6. The mass contours for M^±± = 90GeV at the intermediate scale and fi = 0.25 are shown. The fi mass contours 
(dashed line) from left to right depict 45, 60, 80 and 100 GeV masses, the A ±:t mass contour (solid lines) depicts 88 GeV, the 
in mass contours (dot dashed line) from left to right depict 100 and 140 GeV masses. 



The A will decay into a pair of rs (since the coupling to the other leptons are suppressed). The r's have high px- 
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We have 4 ts in the final state. The like charged ts originate from the same vertex. 

The SM background for this process can come from the pair production of and the subsequent decay of each 
Z a to a t + t - pair (in all these cases oppositely charged rs come from the same vertex). This event will have no 
missing energy, but the rate is small due to the small branching ratio for Z° — > t + t~ (a.B 2 ~ 10~ 3 pb). The 4t 
background can also come from the production of Zj*, 7* converting into t + t~ and the Z decaying into t + t~ pair. 
Another source can be the production of two virtual photons along with e + e~ and each photon converting into t + t~. 
However, in both the above processes, the cross section is ~ 10~ 3 pb. Thus, the SM background is negligible small for 
the 4r signal. 

The A will primarily decay into f\ and r (almost 100%). The fi will then decay into r and Xi (missing energy) 
(100%). The final state has 4 r's plus missing energy. Two of the rs have high pt and these originate from the decay 
of lighter stau. This kind of 4 r plus missing energy signal also originate from the Xi P a i r production in the GMSB 
scenarios where lighter stau is the NLSP. But there is a subtle difference in the final state which we will discuss in 
the next section. 

The chargino pair production can also give rise to 2r plus missing energy states. Since the staus are much lighter 
than the other sleptons, the chargino will primarily decay into fi and v T in the leptonic decay channel. On the other 
hand, the chargino decays into e's and /z's in MSSM. The charginos can be pair produced at LEP II and Tevatron. The 
production of chargino and the second lightest neutralino (this crosssection is larger than the chargino pair production 
) at the Tevatron will also give rise to lots of high pr taus in the final state. The second lightest neutralino (x§) 
primarily decays into t and fj and T\ then decays into a tau and the lightest neutralino. Altogether, there can be 3ts 
along with missing energy. Since X2 m &ss is much larger than the f\ mass, all the three rs will have high px- 

The A bosons and the A-inos can be pair produced at the LEPII and at the Tevatron as well. The productions of 
A's (fermion and boson) at LEP II and A scalar at Tevatron have been considered in the references |l7| ]. We give the 
production cross-sections of the A bosons at the LEP II and at the Tevatron in Figures 7 and 8. The A-ino is pair 
produced at the LEP II and at the Tevatron via Z, 7 exchange. Usually there is also a selectron mediated t-channel 
contribution in the case of A-ino production at the e + e~ collider. However the contribution from this diagram in this 
model is negligible since the A coupling to the first generation leptons is very small. The production crosssections 
are larger than the scalar counterparts for the same mass. In Figure 9 we show the production of A at LEPII for 
yfs = 182 GeV and 194 GeV. We can see that the A-ino mass of 95 GeV gives rise to a cross section of 3 pb at LEPII. 
Thus, if the A-inos are produced at LEPII the cross section will be quite large. 

In Figure 10 we show the production cross-section at the Tevatron for y/s = 2TeV. The production crosssection 
is about 0.7 pb at y/s = 2TeV for A-ino mass of 95 GeV and the cross section is about 0.6 pb at y/s — 1.8TeV for 
the same A-ino mass. Hence with llOpfe -1 of already accumulated luminosity, the number of events are ~ 66. Since 
the final state are pure r leptons, detection is difficult. But we will have some events left even after taking the r 
detection efficiency to be small. We urge the experimentalists to look for the r signals in the data which has been 
already accumulated and also in the data that will be generated in the future runs. 
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FIG. 7. The A ±± pair production GeV at the cross section at LEPII. The solid and the dotted line correspond to the cross 
sections at the center of mass energies to be 194 GeV and 182 GeV. 
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FIG. 8. The A ±=t pair production cross section at Tevatron (center of mass energy is 2 TeV) 
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FIG. 9. The A ±=t pair production cross section at LEP II. The solid line and the dotted line correspond to the cross section 
at the center of mass energies to be 194 GeV and 182 GeV respectively. 
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FIG. 10. The A pair production cross section at Tevatron (center of mass energy is 2 TeV) 

V. GAUGE MEDIATED SUSY BREAKING SCENARIO 

In the previous section we assumed that the scale at which supersymmetry is broken is higher than the Wr scale. 
However this need not be the case and in particular there has recently been a lot of interest in theories where gauge 
interactions are the mediators of supersymmetry breaking at a relatively low scale Q . Let us discuss the implications 
of this scenario for our model. The soft SUSY breaking terms are now generated explicitly only at the scale at which 
the messenger fields are integrated out, and are not explicitly present at the Wr scale. Since they are generated by 
loop graphs involving the gauge bosons of the residual symmetries, their form will be such as to respect only the 
surviving gauge symmetries. We will show that this difference has consequences for phenomenology. Let us assume 
for simplicity that the messenger sector consists of a vectorlike isosinglet pair of fields (charge -1/3) Q © Q and a 
vectorlike weak isodoublet pair L L. 

As shown in 0, in this case, the scalar A mass has two contributions at the SUSY breaking scale of 100 TeV or 
so: one coming from the two loop gauge contributions in the usual manner (for a review see JTi|]) and another coming 
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from the higher dimensional terms. On the other hand, the A mass gets contribution only from the latter kind of 
terms. We therefore expect that in the gauge mediated scenario, the A will be lighter of the two particles. 
The soft SUSY breaking gaugino and the scalar masses at the messenger scale M are given by 

and 

m 2 (Af) = 2n / (A) ^ h Q (^'a 2 

where oti (i = 1 — 3) are the three SM gauge couplings and ki — 1, 1,3/5 for SU(3), SU(2), and U(l), respectively. 
The Ci are zero for gauge singlets, and 4/3, 3/4, and (Y/2) 2 for the fundamental representations of SU(3) and SU(2) 
and U(l)y respectively (with Y defined by Q = 1% + Y/2). Here n corresponds to n((Q, L) + (Q, L)). g(x) and f(x) 
are messenger scale threshold functions with x = A/M. 

We calculate the SUSY mass spectrum using the appropriate RGE equations Jll| with the boundary conditions 
given by the equations above and we vary five free parameters A, M/A, tan/3, n and the sign of fj, (fj, is the coefficient of 
the bilinear Higgs term in the superpotential). We first run the Yukawa couplings (along with the three new couplings 
/i,2,3) and gauge couplings from the weak scale upto the GMSB scale. At the GMSB scale we use the boundary 
conditions and then use the necessary RGEs for the soft SUSY breaking masses in order run down to the weak scale. 

The CLEO constraint on the b — > 57 rate restricts /i < fliSfl . In the absence of late inflation, cosmological 
constraints put an upper bound on the gravitino mass of about 10 4 eV H, which restricts M/A = 1.1 - 10 4 . In 
the figures we show the results for n=l, but we discuss the other values of n also. For reasons discussed before, we 
assume that /i^ are small ( ~ 0.05), but /3 is larger. We show our result for ~ 0.5 and /a ~ 0.25. We also vary 
M A between 70-120 GeV at the GMSB scale. 

The gravitino is always the LSP. In the usual GMSB case the Xi an( i the fi fight for the NLSP spot. In this 
model the A ±=t also joins the race to become NLSP. The third generation lighter stau mass gets affected due to the 
presence of the additional large coupling fa. Thus the f\ is much smaller compared to the conventional GMSB case 
for the same parameter space. Consequently the lighter stau will be the NLSP for a wider region of parameter space 
compared to the lighter neutralino. 

In figure 11a we have shown the mass contours of fi (solid line), Br (dotted line), \i (dash-dotted line) and the 
chargino masses (dashed line) for M = 1.1 A, n=l and M A = 90 GeV at the GMSB scale (94 GeV at the weak scale). 
We also show the contour along which the lighter stau mass and the neutralino mass are same (thick solid line). The 
region above the contour has lighter stau as the NLSP. We see that only a small region for tan/3 3-15 and A 40-60 
TeV has Xi as NLSP. When Xi is the NLSP, it decays into a photon and a gravitino. If Xi is P&i r produced at LEPII, 
the final state has 77 plus prp . The photons are hard and easy to detect. Already we have bound on the Xi mass of 
around 80 GeV at LEP II |2(jj provided the selectrons are not too heavy. We can see that if we use Xi mass bound as 
80 GeV, the region left out in fig 11a where Xi is the NLSP is very small. The region where fi is the NLSP, which is 
the dominant region, stau decays into a high px r and a gravitino (missing energy) . So far there is not much bound in 
these regions, other than the f\ has to be larger than 57 GeV. In the figures we have shown the stau mass contours of 
45, 60, 80 and 100 GeV. The r\ mass contours have large dependence on the tan/3, the t\ decreases with the increase 
in tan/3. The Br, chargino and the Xi mass contours do not have much tan/3 dependence. 

In figure lib we show the same mass contours in the plane of A and M/A for tan/3 = 10. The contours along 
which fi mass is equal to the Xi mass form an envelope. Within the envelope Xi is NLSP. The &r masses increase as 
the M/A ratio increase. The couplings /1.2 are small to affect the Sr or the jiR masses (selectron and smuon masses 
are almost same). The lighter stau mass however decreases with the increase in the ratio of M/A (subtractive effect 
from the soft scalar masses due to the presence of the new coupling overcome effects coming from the gaugino masses 
in the RGE). In the conventional GMSB model the stau mass increases with the increase in the M/A ratio. With 
the improvement of the f\ mass bound much more parameter space will be ruled out in the higherM/A ratio. The 
Xi and the chargino masses initially decrease with the increase in the ratio of M/A due to the threshold corrections. 
The 80 GeV Xi bound rules out much of the parameter space in the envelope wherex? is the NLSP. The M A±± varies 
between 94-108 GeV at the weak scale. Unlike the SUGRA scenario, the A ±± mass is much larger than the fermionic 
part at the weak scale and its mass, in the full range shown in the figure, varies from 170 -570 GeV (due to the large 
soft SUSY breaking contribution at the GMSB scale). Thus, the allowed parameter space in GMSB scenarios are 
much more than the SUGRA scenarios. 

In figures 12 and 13, we have shown the mass contours in the A — M/A plane for M A±± = 70 and 110 GeV at 
the GMSB scale. The tan/3 is chosen to be 10. When M A±± = 70 GeV, the envelope is larger, since f\ mass has 
less subtractive contribution from the delta mass. On the other hand when M A±± = 110 GeV, the envelope shrinks. 
The envelope can also increase in the size if we have smaller coupling . We show the effect of a smaller fa in figure 
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14. This figure looks more like the conventional GMSB model. We see that there is no parameter space where f\ is 
NLSP. But if we increase tan/3 we will hit the region where f\ is the NLSP. The f\ mass in this case, as expected, 
increases with the increases in the ratio of M/A. So far in all these figures we have used n=l. The effect of larger 
values of n can easily be surmised from the mass formula. As n increases gaugino masses increase proportionally, on 
the other hand the scalar masses increase as \fn. Hence stau mass is the NLSP for a even wider region of parameter 
space. The lighter selectron's mass becomes closer to the lighter stau mass and lower than the neutralino mass. 

Let us now discuss the signals. At LEP II, the main production processes are the Xi pair, the A ±=1= pair, the en 
pair and the f% pair. 

In the case of A ±± pair production, each A ±=t wll decay primarily into a fx and a tau (both having same sign of 
charge). The other decay modes involving e.g. a electron and a en and a muon and a smuon are suppressed primarily 
because of small /1.2 couplings. In the parameter space when stau is the NLSP, A ±± decays into a fi and a r and 
fi — > tG. Thus, the final state has 4r plus missing energy and out of the 4 rs, two have high px (higher px than the 
SUGRA case). These two high px rs have opposite sign electric charges. 

The Xi, Pair production also gives rise to 4 rs plus missing energy (since each Xi —> r i r ), with two of these rs having 
high pt [pi| . However there is an essential difference between the signal in this case and in the previous case. Since 
the neutralino is a majorana particle, the two high px rs can have same or opposite sign (with equal probability) for 
the electric charges and thus providing a way to discrminate between the two cases. 

In the case of en pair productions, each selectron will either decay into a electron and a x? or a A and a positron 
(if the selectron mass is higher than the Xi mass and the A ±=t mass). Both Xi and A ±± decay to 2r plus missing 
energy. The final states in either case will have 2e4r plus missing energy with two of the rs having high px- Their 
relative sign will determine the decay channel of the selectron. If, however, the en mass is lower than both the Xi mass 
(large n case) and the A ±=t mass, the selectron can decay into an electron and a gravitino or via offshell production 
of the Xi or A ±:t . The Xi and A ±=t then convert into 2r plus missing energy. In the former case the final state of 
the selectron pair production has 2e plus missing energy and in the later case the signal is 2e4r plus missing enrgy. 
Depending on the parameter space, these offshell decay modes can be comparable or greater than the onshell decay 
mode p2[ . In the case of smuons the electrons in the final states will be replaced by the muons. 

When Xi is the NLSP, the A ±:t will decay into a fi and a r (100%). The fi will decay into a Xi and a r. The 
Xi then decays into a photon and a gravitino. All the branching ratios are 100 %. The final state from the A ±=t 
pair production will have 274r plus missing energy-which is a spectacular signal and very hard to miss. This signal 
will appear along with the electrons or the muons in the case of selectron or smuon production, where the selectron 
decays as described above i.e. en — > eA ±± . 

In the conventional GMSB models when Xi is the NLSP case, one gets at most 2 leptons along with 2 7 (through 
the slepton productions whose cross section is much smaller than the delta-ino production cross section) at the LEP 
II and at the Tevatron or 3 leptons plus 2 photons (chargino-second lightest neutralino production) at the Tevatron. 
Thus the signal 274rplus missing energy will clearly distinguish this model from the ordinary GMSB models in the 
parameter space where Xi is the NLSP. It may also happen that the A ±± mass is smaller than the f\ or the Xi mass. 
In that case the A ±=t will decay into a r and a virtual fi which will convert into a r plus gravitino or if the Xi mass 
is lower than the fi mass, then the t\ will decay into a r and a Xi- The Xi will then convert into a photon and the 
gravitino. Thus the final states are same as discussed in the cases when Xi is the NLSP or fi is the NLSP. 
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FIG. 11. a) The mass contours for M^±± = 90GeV at the GMSB scale and fz = 0.5 are shown. The fi mass contours 
(solid line) from left to right depict 45, 60, 80 and 100 GeV masses, the Er mass contours (dotted lines) depict 80 (along the 
left axis), 95 and 115 GeV masses, the Xi (dot-dashed) mass contours depict 60, 80 and 100 GeV masses, the lightest chargino 
(x ) mass contours (dashed) depict 94, 100 and 130 GeV masses b) The mass contours for M^±± — 90GeV at the GMSB 
scale and fz = 0.5 are shown. The fx mass contours (solid line) from left to right depict 45, 60, 80 and 100 GeV masses, the 
in mass contours (dotted lines) depict 80, 95 and 115 GeV masses, the X? (dot-dashed) mass contours depict 60, 80 and 100 
GeV masses, the lightest chargino mass contours (dashed) depict 88, 100 and 130 GeV masses. The thick solid lines in both 
the figures depict the contour along which the lighter stau mass equals to the lightest neutralino mass. 
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FIG. 12. The mass contours for M^±± = 70GeV at the GMSB scale and fz — 0.5 are shown. The f\ mass contours 
(solid line) from left to right depict 45, 60, 80 and 100 GeV masses, the Br mass contours (dotted lines) depict 75, 95 and 
115 GeV masses, the x? (dot-dashed) mass contours depict 60, 80 and 100 GeV masses, the lightest chargino mass contours 
(dashed) depict 70, 100 and 130 GeV masses. The thick solid lines depict the contour along which the lighter stau mass equals 
to the lightest neutralino mass. 
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FIG. 13. The mass contours for M^±± — HOGeV at the GMSB scale and f% — 0.5 are shown. The ti mass contours (solid 
line) from left to right depict 45, 60, 80 and 100 GeV masses, the en mass contours (dotted lines) depict 95 and 115 GeV masses, 
the Xi (dot-dashed) mass contours depict 60, 80 and 100 GeV masses, the lightest chargino mass contours( dashed)depict 100 
and 130 GeV masses. The thick solid lines depict the contour along which the lighter stau mass equals to the lightest neutralino 
mass. 
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FIG. 14. The mass contours for M^±± — 90GeV at the GMSB scale and fa = 0.25 are shown. The ?i mass contours (solid 
line) from left to right are 80 and 100 GeV, the en mass contours (dotted lines) are 75, 95 and 115 GeV, the x° (dot-dashed) 
mass contours are 60, 80 and 100 GeV, the lightest chargino mass contours( dashed)are 70, 100 and 130 GeV. 



VI. CONCLUSION 

To summarize, wc have studied the phenomenological implications and the collider signatures of the remnant doubly 
charged Higgs boson A ±:t and its fermionic partner the Deltaino A ±± . We find that existing limits on lepton flavor 
violation and muonium-anti-muonium transition can be used to conclude that the dominant coupling of these particles 
are diagonal and mostly to the third generation of leptons. They of course are forbidden from coupling to the quarks 
by electric charge conservation. The effect of the dominant third generation lepton coupling has the consequence that 
the fi mass is much smaller than the other slepton masses. This gives rise to mulitiple r enriched signal at LEP II 
and Tevatron. The fermionic as well as the bosonic partner of the doubly charged Higgs bosons can be produced at 
the present colliders and the signals contain 4 r (plus 2j in some scenarios of the GMSB version of the models) with 
and without missing energy. This signal is found to be detectable for reasonable range of mass values of the particles 
and could be used to test the supersymmetric left-right models of the type discussed here or to restrict the allowed 
parameter range of the model. Most of the allowed parameter space can be searched in the existing colliders. We 
therefore urge the experimentalists to analyze the tau events in the existing as well as in the future data. 
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VII. APPENDIX 



The standard MSSM RGEs for the Yukawa couplings and the soft SUSY breaking terms will be modified due to 
the presence of the new couplings and the fields,. There will be new RGEs for the additional fields and the couplings. 
The new interaction involves the third generaion righthanded leptons and the new field A ±=t . We will keep only the 
/ 3 coupling in the RGE's (in our numerical calculation we use all of them). The new RGEs and the modified ones 
are listed below: 

2VX T = A T (- environ) + X 2 b + 4A 2 + 4/ 3 2 ), (9) 



where Cr^ = 9/5, 3, and T> 



where Cf 3i = 36/5,0,0. 



IOtt- 



2 dt ' 



2Vf 3 = / 3 (- J2 Cfaii^m) + 10/ 3 2 + 4A 2 ), (10) 



r>m 2 c = -— irayM 2 + 2A 2 (m 2 c + m 2 + M% x + A 2 T ) (11) 



+ 4/ 3 2 (m 2 , ±± +M 2 c +V), 
192 

Vm 2 A ±± = — iraxMf + 2f 3 2 {m 2 A±± + 2m 2 .. + A 3 2 ), (12) 

o 

We do not include the effects of the couplings to the other generation f\ ; % since they are very small. 

VA T = {J2 CniwaiMi + X 2 A b + AX 2 T A T + 4/ 3 2 A 3 ), (13) 



VA 3 = ((47rai)y M! + 10/ 3 2 A 3 + 4A 2 ^ T ), (14) 



2VM K±± = M A±± (-(47rai)— My + 4/ 3 2 ). (15) 

o 
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